RNA polymerase motor on DNA track: 
effects of interactions, external force and torque 
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RNA polymerase (RNAP) is a mobile molecular workshop that polymerizes a RNA molecule by 
adding monomeric subunits one by one, while moving step by step on the DNA template itself. 
Here we develop a theoretical model by incorporating the steric interactions of the RNAPs and 
their mechanochemical cycles which explicitly captures the cyclical shape changes of each motor. 
Using this model, we explain not only the dependence of the average velocity of a RNAP on the 
externally applied load force, but also predict a nonmotonic variation of the average velocity on 
external torque. We also show the effect of steric interactions of the motors on the total rate of 
RNA synthesis. In principle, our predictions can be tested by carrying out m-vitro experiments 
which we suggest here. 

PACS numbers: 87.16.Ac 89.20.-a 
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I. INTRODUCTION 

Molecular motors m in living cells arc either proteins 
or macromolecular complexes made of proteins and ri- 
bonucleic acids (RNAs). Like their macroscopic coun- 
terparts, these motors perform mechanical work, while 
translocating on a filamentous track, by converting in- 
put energy which is often supplied as chemical energy 
HIjIIjHI- In this paper we study a special class of motors, 
called RNA polymerase (RNAP) which play crucial roles 
in gene expression [J. 

Transcription of a gene encoded in the sequence of nu- 
cleotides in a specific segment of a DNA is carried out b 
RNAP motors which treat the DNA as a template 
An RNAP is more like a mobile workshop that performs 
three functions simultaneously: (i) it decodes the genetic 
message encoded in the template DNA and selects the 
appropriate nucleotide, the monomeric subunit of RNA, 
as dictated by the template, (ii) it catalyzes the addition 
of the monomeric subunit thus selected to the growing 
RNA molecule, (iii) it steps forward by one nucleotide on 
its template without completely destabilizing the ternary 
complex consisting of the polymerase, the template DNA 
and the product RNA. The free energy released by the 
polymerization of the RNA molecule serves as the input 
energy for the driving the mechanical movements of the 
RNAP. Therefore, these enzymes are also regarded as 
molecular motors 

During the transcription of a gene, the collective move- 
ment of the RNAPs on the same DNA track is often re- 
ferred to as RNAP traffic because of its superficial sim- 
ilarity with vehicular traffic [1, Q. The beginning and 
the end of the specific sequence corresponding to a gene 
are the analogs of the on-ramp and off-ramp of vehicu- 
lar traffic on highways. The average number of RNAPs, 



which complete the synthesis of a RNA molecule per unit 
time interval can be identified as the flux in RNAP traf- 
fic. Note that flux is the product of the number den- 
sity and average velocity of the motors. Thus, fiux in 
RNAP traffic is identical to the average rate of synthe- 
sis of the corresponding RNA. Using the terminology of 
traffic science Q, we'll call the relation between the flux 
and the number density of the motors as the fundamen- 
tal diagram. The fundamental diagram is an important 
quantitative characteristic of traffic flow. 

The dependence of the velocity of the motor on an 
externally imposed load (opposing) force is called the 
force- velocity relation which is one of the most important 
characteristics of a molecular motor. The force-velocity 
relation for RNAP motors have been measured by carry- 
ing out single molecule experiments El • However, to 
our knowledge, the response of an RNAP motor to an ex- 
ternally applied torque has not been investigated so far. 
The effects of steric interactions of the RNAP motors on 
their dynamics has been studied only in a few experi- 
ments H [ll El, [H, [3 ; but, none of these addressed 
the question of the nature of the overall spatio-temporal 
organization of the RNAP motors in RNAP traffic. 

The traffic-like collective dynamics of cytoskeletal 
molecular motors [13, [H, [ilj20j2l|j2i and that of 
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ribosomes on mRNA tracks 

[sol , [sil . [3^ . [33I [s^ l have been investigated theoretically 
in the physics literature. However, so far, RNAP traf- 
fic has received far less attention [H, [Sy, [s^] In this 
paper we develop a model that captures not only the 
steric interactions between the RNAPs, but also sepa- 
rately the biochemical reactions catalyzed by an RNAP 
and the cyclic shape changes it undergoes during each 
mecliano-chemical cycle. This model may be regarded 
as a "unified" description in the sense that the same 
model describes the single RNAP properties (e.g., the 
force- velocity and torque- velocity relations) as well as the 
collective spatio-temporal organization (and rate of RNA 
synthesis) . 



II. BASIC MECHANO-CHEMISTRY AND THE 
MODEL 

The main stages in the polymerization of polynu- 
cleotides by the polymerase machines are common: 

(a) initiation: Once the polymerase encounters a specific 
sequence on the template that acts as a chemically coded 
start signal, it initiates the synthesis of the product. The 
RNAP, together with the DNA bubble and the growing 
RNA transcript, forms a "transcription elongation com- 
plex" (TEC). This stage is completed when the nascent 
product RNA becomes long enough to stabilize the TEC 
against dissociation from the template. 

(b) elongation: During this stage, the nascent product 
gets elongated by the addition of nucleotides; during 
elongation [sl], each successful addition of a nucleotide 
to the elongating mRNA leads to a forward stepping of 
the RNAP. 

(c) termination: Normally, the process of synthesis is ter- 
minated, and the newly polymerized full length product 
molecule is released, when the polymerase encounters the 
terminator (or, stop) sequence on the template. 

In this paper we are interested mainly in the elongation 
of the mRNA transcripts. 
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FIG. 2: A schematic representation of the mechano-chemical 
cycle of each RNAP in our model in the elongation stage 
is shown in (a) (see the text for explanation). The labels 
— + in (b) denote the nucleotides on the DNA 
template. Thus, the horizontal direction in (b) denote po- 
sitions of the RNAP while the vertical direction, at a given 
position, describe the corresponding "chemical" states. 




FIG. 1: (Color online) A schematic depiction of the essen- 
tial architectural feature of a RNAP that we capture in the 
model proposed in this paper. The RNAP can switch be- 
tween "open" and "closed" shapes cyclically in each mechano- 
chemical cycle of its operation. The extent of "opening" is 
measured by the angle 6 which can be manipulated also by 
externally applied torque r. 

One of the key experimental observations is that an 
RNAP cycles between a "closed" and an "open" shape 
during each mechano-chemical cycle (see figlHand figH^) 
[sol , lio . l4ll . . Recognition of the correct incoming nu- 
cleotide (dictated by the template) and its binding with 
the catalytic site on the RNAP leads to its closing. It re- 
mains closed as long as the hydrolysis of the NTP takes 
place. The nascent RNA is elongated by one monomer 
on completion of this reaction which also produces py- 
rophosphate. Then, the shape of the RNAP switches 
back to its open shape from the closed one thereby facil- 



itating release of the pyrophosphate. In the open shape, 
the RNAP is weakly bound to its DNA track and, in 
principle, can execute Brownian motion. 

To our knowledge, none of the models of RNAP traf- 
fic reported earlier, explicitly capture these cyclic change 
of shape of the polymerase motor. Liverpool ct al.[3H| as 
well as Klumpp and Hwa [IJ treat each RNAP as a rigid 
rod and describe its translocation from one nucleotide to 
the next on its template in terms of a single parame- 
ter. Thus, both these models capture the effects of the 
mechano-chemistry of indivisual RNAPs by an effective 
"hopping" rate contstant. In an attempt to capture the 
most essential aspects of the mechano-chemical cycle of 
individual RNAP motors, we have assigned two possible 
"internal" states to an RNAP at each spatial position on 
its track [s^. However, each of these two states is, in 
reality, not a single bio-chemical state; a sequence of bio- 
chemical states which interconvert sufficiently rapidly are 
collectively represented by one "internal" state whereas 
the transition from one internal state to the other is the 
slowest of all the transitions and, therefore, rate-limiting. 
This minimal model neither distinguishes between purely 
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chemical transitions pure shape changes, nor does it as- 
sign distinct internal states to the "open" and "closed" 
shapes of the RNAP motor. Therefore, in this paper, we 
extend our earlier 2-state model to a 4-statc model. 

We shall use the terms RNAP and TEC interchange- 
ably. We extend the 2-state model of ref. jsH] by assigning 
four possible distinct "internal" state to the RNAP to 
capture its transitions between open and closed shapes. 
The composition of these four states and the spatial po- 
sition of the RNAP as well as the transitions between 
these states are shown in figl^l^b). 

The subscript n of the symbol "RNAP" denotes the 
position of the RNAP motor on its track, measured 
from a nucleotide that marks the starting point of 
transcription by the RNAP. The two superscripts (o) 
and (c) of "RNAP" refer to its "open" and "closed" 
shapes, restectively. The four distinct states in each 
cycle are labelled by the integers 1,2,3,4. Let us begin 
with the state 1 , where the RNAP is open and is located 
at n so that the length of the corresponding elongating 
mRNA transcript is also n. Arrival of the correct NTP 
subunit causes the transition to the state 2 where the 
open RNAP, now located at the new position n -I- 1, 
forms a complex with the mRNA transcript of length n 
and the newly arrived NTP subunit. Next, the shape of 
the RNAP changes from open to closed form and the 
new state of the system is labelled by the integer index 
3. Then, the closed RNAP catalyzes the formation of 
a new covalcnt bond between the mRNA transcript 
and the newly arrived NTP subunit which leads to the 
elongation of the transcript from n to n -\- \ and the 
corresponding transition from state 3 to state 4. Finally, 
the RNAP opens again and releases the pyrophosphate 
{PPi)] this transition takes the open RNAP to the state 
1 at the position n + 1. The 4-state model reduces to our 
earlier 2-state model [s^ in the limit ^23 00,(^34 00 



III. RNAP TRAFFIC UNDER PERIODIC 
BOUNDARY CONDITION 

In our model a one dimensional lattice of L sites rep- 
resents the DNA template. Each site of these lattice 
corresponds to a single base pair. The total number of 
RNAP motors which move simultaneously on this tem- 
plate (and transcribe the same gene sequentially) is N; 
the linear size of each RNAP is r in the units of base 
pairs. Therefore, the number density of the RNAP mo- 
tors is p = N/L whereas the corresponding coverage den- 
sity is pcov = Nr/L. The instantaneous spatial position 
of a RNAP is denoted by the leftmost site of the lattice 
it covers at that instant of time; however, it also "cov- 
ers" all the adjacent r — 1 on the right side of of the site 
which denotes its position. In order to capture the steric 
interaction (mutual exclusion) of the RNAPs, no site of 
the lattice is allowed to be be covered simultaneously by 
more than one RNAP. Thus, this model of RNAP traffic 
may be regarded as a totally asymmetric simple exclusion 
process (TASEP) [i^, |3| for hard rods each of which can 
exist in one of its four possible "internal" states at any 
arbitrary position on the lattice. 

Let P^(i,t) denote the probability that there is a 
RNAP at the spatial position i and in the chemical state 
/i at time t. At any arbitrary site n and time these 
probabilities must satisfy the normalization condition 
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P(n,i) = 5]P,,(n,t) = -^-p (1) 

^^=l 

Under the mean-field approximation, the master equa- 
tions for P^{i,t) are as follows: 



dPi {n, t) 
dt 



dP2{n,t) 
dt 



= Pi(n + l,t) tjJi Q(n + 1 - r\ n + l ) + P2{n + l,t) Q(n + 1 - r\ n + l ) 
-f P4{n,t) uj^i — Pi{n,t) uj\iQ{n — r\n) ^- u}{2 Q{n\n + r) 



Pi{n - l,t) Lui^ Q( n-l \n - 1 + r) - P2{n,t) [ uj^^i Q{n ~ r\n) + W23] 



dP3{n,t) ^ 
dt 

dP4,{n,t) 
di 



P2{n,t) UJ23 - P3in,t) W34 



= Pz{n,t) W34 - P4 W41 



(2) 
(3) 
(4) 
(5) 



where Q{i\j) is the conditional probability that, given an RNAP in site i, site j is empty. We calculate Q{n\n + r) 
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following the same steps as we did in ref. [3^ , getting 

I — pr 



Q{n\n + r) = Q{n\n+_r) 



1 



p- pr 



(6) 



In the steady state, all the P^(n,t)'s become indepen- 
dent of time. Moreover, because of the periodic boundary 
conditions, all the sites are equivalent so that the steady- 
state probabilities P^(ri,)'s arc also independent of site 
index n. Solving Eqn. (|2l3l4l5p . together with normal- 
ization condition ((T|), in the steady state under periodic 
boundary conditions, we get 



Pi 
P2 
P3 
Pa 



U>il + a;34 0141 ^21*3 fry\ 

p 1? (7) 
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W34 UJil Lol^Q 
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mi iOi2 ^23 Q 
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(^(2 ^23 ^34 Q 
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(8) 
(9) 
(10) 



Where; 



K = U!23i^340J41 + 0J34(^4li^2lQ + ^3i^il^l2Q 

+ UJAiUj{^UJ23Q + ^i2^23^3AQ (H) 

The corresponding flux J is given by 



J 



Q 



K 



^3i ^41 



^12 ^23 — ^\l '^23 ^ t^ll ^21 Q 



(12) 



Hence, the average velocity y of a single RNAP is given 
by 



y = J/p W34 a;4i 



W12 W23 - a;23 - uj\-^ UJ21 



For numerical calculations, we use the rate constants 
extracted by Wang et al.[4^ from the empirical data. 
In the absence of external force and torque, the typical 
numerical values of these rate constants are as follows: 



The average velocity V depends on the number density 
p through the p-dependence of Q. 

In the following subsections we'll use the formula 
in the regime of sufficiently low coverage density pcov 
to predict the dependence of V of single RNAP motors 
on external force and torque. Such dependences of V 
on external force and torque can be probed by carry- 
ing out single-molecule experiments in-vitro. In order 
to predict the total rate of RNA synthesis at coverage 
densities where steric interactions between the RNAPs is 
not negligibly small, we'll use the expression Thus, 
our model may be regarded as a "unified" description of 
transcription in the sense that it can account for prop- 
erties of single RNAP motors as well as their collective 
behaviour. 
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• [NTP], with cjf" = 10*^ M^^ ■ s-^ 



^23 = 100.0 s"^ 
tJ34 = 1000.0 fi-i 
CJ41 = 31.4 

(14) 

So far as the two new rate constants ^23 and 0^41 are 
concerned, we investigate their effects on the rate of tran- 
scription by varying their magnitudes over a wide range. 



A. Effects of externally applied force and torque 

In this section we examine the effects of external force 
F and torque r on the rate of RNA polymerization. From 
the force- velocity relation, we extract the stall force Fs, 
the load force which just stalls the a single RNAP motor. 
In order to extract these single RNAP properties from 
our model, we first use the formula \Vi\ in the regime of 
extremely low coverage density of the RNAP motors. 
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FIG. 3: The average speed of a single RNAP motor is plot- 
ted against the load force (a) for three different values of the 
NTP concentration. The lines have been obtained from the 
analytical expression (|13|l in extremely low number density of 
the RNAP motors. The discrete data points have been ob- 
tained by computer simulations of the model under identical 
conditions. 



1. Ejfects of external force 

Following the standard practice in the master equa- 
tion approach to molecular motors 0] for capturing the 
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FIG. 4: The average speed of a single RNAP motor is plotted 
against the torque (a) for three different values of the NTP 
concentration for fixed CJ23 = 100s~^, (b) for four different val- 
ues of the rate constant LJ23 , keeping NTP concentration fixed 
at 10~^M. The lines have been obtained from the analyti- 
cal expression (|13p at extremely low coverage density of the 
RNAP motors. The discrete data points have been obtained 
by computer simulations of the same model under identical 
conditions. 



effects of external force on the mechano-chemistry [46|, 
we multiply the "mechanical" steps of the cycle by an 
appropriate Boltzmann-like factor. More precisely, we 
assume that the load force F significantly affects only 
those steps of the mechano-chemical cycle of an RNAP 
which involve its mechanical movement, i.e., forward or 
backward movement in real space. Therefore, we assume 
that force-dependence of these rate constants are given 
by 
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exp{-F S/keT) 
exp{F S/keT) 
exp{F S/keT) 



where 01(2,1021 '^11 the corresponding values of 
the rate constants in the absence of the load force. In 
equation (|15p the symbol 6 ~ 0.34 nm is the typical 
length of a single nucleotide. None of the rate constants 
other than the three listed in (fT5|) are affected by the load 
force F. 

The force-velocity relation of individual RNAP mo- 
tors in our model is shown in figO for three different 
concentrations of the NTPs. The average velocity de- 
creases monotonically with increasing load force; the con- 
vex shape of the load-velocity curves are very similar to 
those reported by Wang et al. [1^ in their single RNAP 
model. In this model, the higher is the NTP concentra- 
tion the larger is the stall force Fs. Moreover, for a given 
load force F < Fs, the average velocity of the RNAP is 
larger at higher NTP concentration. 



2. Effects of external torque 

Now we consider the effects of an externally imposed 
torque which assists opening, but opposes closing of the 
paZm-like shape of an isolated RNAP. These compet- 
ing effects of the same torque on the two steps of each 
mechano-chemical cycle of individual RNAP motors has 
nontrivial effects on the rate of RNA synthesis. The 
above mentioned effects of the torque r on opening and 
closing of RNAP are captured by the following choice of 
the corresponding rate constants 



W23(t) = L02S. exp{-T 6/kBT) 
W41 (r) = LJ41 exp{T d/ksT) 



(16) 



where LU23 and ^41 on the right hand sides of the equa- 
tions are the rate constants in the absence of external 
torque. 

The average velocity of the RNAP motors are plot- 
ted against the torque in fig[4]for different sets of values 
of the model parameters. The most dramatic effect of 
the torque is the nonmonotonic variation of V with t 
in several experimentally accessible regimes of parame- 
ter values. Increasing torque increases the opening rate 
and decreases the closing rate. As long as the open- 
ing rate is still lower (and, hence, rate limiting), the 
velocity increases with increasing torque. A peak ap- 
pears where closing becomes the rate limiting step. In 
order to demonstrate the effect of this crossover from a 
regime dominated by opening to that dominated closing 
of RNAP, we have plotted one curve in figHJ^b) corre- 
sponding to UJ23 = oJii for which the maximum occurs at 
T = 0. 



Effects of steric interactions: flux-density 
relation 



(15) 
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We have plotted mean-field estimate ([T2|) of flux J 
against the coverage density pcov in figEl for three dif- 
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p = 1/2. The physical reason for the asymmetric shape 
of the fundamental diagram in fig|5]is as follows: to get 
maximum current, system should have maximum num- 
ber of particle- hole pairs. For r = 1, maximum num- 
ber of particle hole pairs occur when 1 — p = p, i.e., 
p = p'" = 0.5. For r > 2, the magnitude of p™ can be 
found from the corresponding condition 1 — r = p™, 
i.e., = l/(r -|- 1), which corresponds to the coverage 
density 



r/{r + 1). 



(17) 



FIG. 5: The steady-state flux of the RNAPs, under periodic 
boundary conditions, plotted as a function of the coverage 
density pcov for three difi'erent values of the ratio L023/i^4i- 
The vertical line represents the position of our theoretical 
estimate (|17|l of the maximum flux for fixed value of transition 
rates. 

ferent valus of the ratio 0^23 /w4i at fixed NTP concentra- 
tion of 10~^M. In the traffic science literature @ such 
flux-density relations are usually referred to as the fun- 
damental diagram. The qualitative features of the fun- 
damental diagrams in figO arc similar to those observed 
earlier in the two state model of RNAP traffic [s^ . The 
most notable feature of these fundamental diagrams is 
their asymmetric shape; this is in sharp contrast to the 
symmetry of fundamental diagram of a TASEP about 



IV. RNAP TRAFFIC MODEL UNDER OPEN 
BOUNDARY CONDITION 

For modeling transcription, the open boundary condi- 
tion is more realistic than the periodic boundary condi- 
tions. Under this conditon a RNAP can attach to start 
site(labled as n = 1) with rate LOa provided none of the 
first r sites on the lattice is covered by any other RNAP. 
Similarly, after reaching the site n ~ r, an RNAP leaves 
the DNA track with rate ojp. We calculate the condi- 
tional probability Q same way as we have done for the 2- 
state model of RNAP [s^. In this case, under mean-field 
approximation, the master equations for the probabilities 
Pp(?7, t) are given by 



dPi{n,t) 
Jt 



dP2in,t) 
dt 



dPsjn^t) 
dt 

dPi{n,t) 
dt 



LOo, (l-IZ ^(*)) '5(n-l)+ Lo\i Pi{n + l,t)+ Lol^ P2{n + l,t) 

ujii PA{n,t) - Pi{n,t) H{n - 2) Q- (n) 

wfa Pi{n,t)(^l~ H{n- L)^ <3+W - ujp Pi{n,t) 5{n - L) 

Lo{^ P,{n-l,t) H{n-2) Q+{n-l)- uj^s P2{n,t) 
uj\^ P2{n,t) H{n-2) Q" {n) 
(^23 P2{n,t) -[J34 P'i{n,t) 

W34 P?,{n,t) - LOii Pi{n,t) 

I 



1-Jf(n-L) Q-{n + l 



(18) 

(19) 
(20) 

(21) 



Where H{n) is the Heaviside step function defined by and 



Hin) = 



1 if n > 
if n < 0. 



g-(n) = 



if n < r 
if n > r. 



(22) 



(23) 
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FIG. 6: The steady state current of RNAPs under open 
boundary condition for three different values of ratio UJ23 /con 



We have numerically solved the mean-field equations 
P^ - (PT|) in the steady-state to compute the correspond- 
ing flux of the RNAP motors. These mean-field theoretic 
estimates of flux are compared with the corresponding 
data obtained from direct computer simulations of the 
model. These results are plotted as functions of the rate 
constants Wq. and Wji, respectively, in figs[6] (a) and (b). 
The flux rises with increasing 1^(2, but the rate of rise 
decreases with increasing o^/j- Eventually, the flux satu- 
rates because longer rate-limiting. This trend 

of variation of the total rate of RNA synthesis with the 
concentration of NTP is similar to that observed earlier 
in our 2-state model of RNAP traffic Thus, a 2- 

state model of RNAP would be adequate to capture the 
dependence of the rate of RNA synthesis on the NTP 
concentration. 



SUMMARY AND CONCLUSIONS 



^ some important cyclic shape changes in the mechano- 
^ chemical cycle of each RNAP. The new model predicts 
the effects of these shape changes on the rate of RNA 
^Jl^hesis. Moreover, we have used the same model in the 
extremely low density limit to extract the force-velocity 
relation. Finally, we have demonstrated a novel non- 
monotonic variation of the average speed of the RNAP 
motors (and, hence, a nonmonotonic variation of the rate 
of synthesis of RNA) with the increase of an externally 
imposed torque on the individual motors. Nonmonotonic 
variation of average speed with external load force (not 
external torque) has been observed earlier [13, E^, in 
the case of several DNA polymerases. But, the phys- 
ical origin of this nonmonotonicity, as interpreted in 
ref. [13, S, [4^, is quite different from that responsible 
for the nonmonotonic variation of average speed with ex- 
ternal torque. Moreover, the present version of our model 
does not explicitly capture the effect of possible twisting 
of the DNA by the polymerase. In principle, it should 
be possible to test the new predictions of our model by 
carrying out laboratory experiments in-vitro. 

To our knowledge, it may not be possible, at present, 
to apply a torque that would directly oppose (or assist) 
the opening (or closing) of a RNAP. Therefore, we now 
suggest a possible alternative technique. The nonmono- 
tonic variation of the velocity of the RNAP with external 
torque will result also in the case where the torque op- 
poses, rather than assisting, the opening of the individ- 
ual RNAP motors (i.e., opposite to the direction of the 
torque shown in figH]). Suppose, the finger- like domain 
of each RNAP is genetically cross-linked with the palm 
domain by a protein whose entropic elastic constant is 
significant. When such an RNAP tends to open up, the 
cross-linking protein opposes the opening. By repeat- 
ing the experiment with different cross-linkers, whose en- 
tropic elastic constants are different, one can study the 
variation of the average velocity of the RNAP motors 
with the applied torque. 

In spite of some crucial differences, most of the 
polynucleotide polymerases seem to share a common 
"cupped right hand" architecture [sO]. Therefore, it 
should be possible to make minor modifications in our 
4-state model of RNAP so as to develop similar models 
of other polynucleotide polymerases. 



In this paper we have extended our earlier 2-state 
model of RNAP traffic to a 4-state model so as to capture 
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